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Abstract 

We predict the existence of linear discrete rogue waves. We discuss that Josephson effect is the underlying reason for the formation 
of such waves. We study linear rogue waves in continuous system and present an exact analytical rogue wave solution of the 
Schrodinger-like equation. 


1. Introduction 

Rogue waves, sometimes known as freak waves or extreme 
waves, are doubly localized waves in both space and time and 
appear as a peak on a finite background. They have been rarely 
observed in oceans and but recently realized in a number of 
different physical contexts such as optics |[T], microwaves 
and inhomogeneous media 13. Rogue wave solutions have 
been studied mostly for the focusing nonlinear Schrodinger 
equation (NLS). Three exact analytical solutions of the NLS, 
Peregrine soliton Q, Kuznetsov and Ma soliton 13] |3 and 
Akhmediev breather Q, clearly explains the dynamics of rogue 
waves. In the Peregrine soliton, a slight modulation on a uni¬ 
form background grows until it reaches its maximum value as 
a result of the modulational instability and then decays and 
vanishes. The Peregrine soliton was experimentally realized 
in an optical system in 2010 a and then in a water wave 
tank in 2011 0. The Peregrine soliton is the hrst order ra¬ 
tional solution of the NLS and the second-order rational solu¬ 
tion was studied in ifTOl and observed experimentally in ITTl . 
The ratio of maximum amplitude of the rogue wave to the 
background amplitude is 3 for the first order rational solution 
while it is 5 for the second order one. The extension of rogue 
waves to discrete systems has been investigated by some au¬ 
thors ini [El [H El [II [nl [IS [EIEqIEII . The key idea to 
construct a discrete rogue wave is to use a proper choice of ini¬ 
tial held amplitude. 

As it is noted in a recent review article the originating 
mechanism of rogue waves is still a matter of debate (refer¬ 
ences therein). The modulational instability is the most popular 
one but is restricted to nonlinear systems. The interpretation of 
rogue waves in terms of the two Akhmediev breather collisions 
was also made for the nonlinear system Il23]l . We would like to 
emphasize that rogue waves can occur even in the absence of 
nonlinear interaction. Rogue waves has recently been observed 
in a linear microwave system m- The authors in |l2l discussed 
that random positions of metal cones acts like a ’’bad lens” di- 
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recting waves towards some focal points and the intensity be¬ 
comes hve or more times greater than background level. In a, 
rogue waves in a linear inhomogeneous media were analyzed 
and it was discussed that granularity and spatial inhomogene¬ 
ity are the two joint generators of optical rogue waves. In Il24l . 
linear rogue wave formation in oceans was analyzed. A uni¬ 
formly distributed initial currents is supposed to enter the 2D 
spatial system along the y axis and propagate along x-axis. The 
currents are considered to have a single velocity but a Gaussian 
angular spread ( the components of the initial velocity in the 
X and y directions change with y). It was shown that the cur¬ 
rents focus along the propagation direction and consequently 
the density becomes high enough to be interpreted as a rogue 
event. The focusing effect was also analyzed in ESll . 

Motivated by the investigations of rogue waves in linear sys¬ 
tems and by the question of underlying mechanism for the ap¬ 
pearance of large amplitude waves, we inquire whether rogue 
waves occur in linear domain for both discrete and continuous 
system. In this paper, we show the existence of linear discrete 
rogue waves and present an exact analytical rogue wave solu¬ 
tion for the Schrodinger-like equation. 

2. Linear Discrete Rogue Waves 

The propagation of an optical held in a tight binding waveg¬ 
uide array or time evolution of matter wave in a tight bind¬ 
ing optical lattice can be described by the discrete nonlinear 
Schrodinger equation (DNLS). It is a nonintegrable equation 
and posseses two conserved quantities, the total number of par¬ 
ticles (or power) and the total energy. We consider DNLS to 
study formation and evolution of discrete rogue wave. 

= -7(T'„+, + T'„.i) + gml'¥„ (1) 

where 'P„ is the complex wave function at the n-th waveguide, 
t is time (f is replaced by the propagation direction z in an ar¬ 
ray of waveguide), J is the coupling coefficient between «-th 
waveguide and adjacent waveguides, g is the nonlinear interac¬ 
tion strength. Our aim is to find discrete rogue wave solution 
of the DNLS not triggered by modulational instability. We will 
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firstly consider the non-interacting limit g = 0 and then investi¬ 
gate the effect of nonlinear interaction. 

Generally speaking, a discrete rogue wave is the strong local¬ 
ization over a few lattice sites. In ll26l . it was discussed that 
initial conditions must be chosen properly to observe discrete 
rogue waves. The subsequent dynamics is strongly sensitive 
to initial wave function. We look for an initial wave function 
that leads to rogue waves even in the absence of nonlinear in¬ 
teraction. We propose to start with the following uniform wave 
function 

CO 

'T„(f = 0) = A ^ c,(n - «,)' (2) 

i=l 

where the background amplitude A is supposed to be a real 
number, the coefhcients c, and the shifts n, are constant and 
the initial phase <!>(«) is position dependent. The initial density, 
|T'(f = 0)1^, is uniform along the lattice and this wave function is 
therefore a good candidate to look for rogue waves that appear 
from nowhere. The coefhcients c, are undetermined. There are 
inhnitely many terms in the expression of the phase 0(n). For¬ 
tunately, principal aspects of linear rogue waves are captured 
by assuming that only two of them are non-vanishing: ci 0, 
C 2 ^ 0 and c, = 0 for i > 2. As a special case, the wave function 
becomes an exact solution when C 2 - 0 and ci 0. This solu¬ 
tion is known as modulationially unstable. If C 2 + 0, the time 
evolution changes signihcantly as we shall see below. No exact 
solution is available even in the absence of nonlinear interac¬ 
tion, we therefore perform numerical computation to hnd the 
time evolution of the above initial wave function. In our numer¬ 
ical computation, we solve the equation (1) for the periodical 
boundary condition, 'TAr+i(f) = 'Ti(f), where N is the total num¬ 
ber of lattice sites. Suppose hrst the noninteracting limit, g - 0 
and take 7 = 1, A = 0.4 and N = 200. The Fig-1 plots the time 
evolution of the above initial wave for three different parame¬ 
ters. These hgures clearly show that time evolution is sensitive 
to the exact numerical values of ci and C 2 . The Fig-1.a, Fig-l.b 
and Fig-l.c are for ci - n, C 2 - 0.005 (single peak), ci = 0, 
C 2 - 0.03 (W-shaped) and ci = 0, C 2 = 0.06 (zigzag-shaped), 
respectively. In all the of figures, we see some peaks that show 
strong localizations of particles. More precisely, particles are 
focused towards some focal points and then defocused like light 
in an imperfect lens in optics. The heights of the peaks are 
nearly three times larger than the background amplitude. The 
peaks appear from nowhere and disappear shortly as a result of 
the tunneling of the particles to adjacent lattice sites. Therefore, 
they can be interpreted as discrete rogue waves. They aren’t 
due to modulational instability since the nonlinear interaction 
is assumed to be zero. Let us analyze the three hgures in more 
detail. In the Fig-l.a, we see a single rogue wave. The ratio 
of the absolute of the maximum amplitude to the background 
amplitude is equal to 3. The particles are focused towards the 
center of the lattice. A single peak is formed at around f = 50 
and then decays in a short time. The decaying mechanism is 
in fact defocusing of focused particles due to the tunneling to 
adjacent sites. When the particles are defocused to all sites, the 
system starts to behave chaotically. (In the Fig.l, the chaotic 
behaviors are not shown) Let us now discuss the Fig. Lb and 


Fig.l.c. There are now more than one focal point and conse¬ 
quently more than one peak. Since the particles are distributed 
in the peaks, the heights of the peaks decrease and consequently 
the widths of the peaks increase. The latter one shows us that 
the lifetime of the rogue waves increases. We would like to note 
that rogue waves are formed earlier in the last two figures. As 
in the Fig. La, the particles in each peaks are spread mainly into 
two lobes and the system behaves chaotically at large times. 
To this end, we note that changing the background density A 
doesn’t change the rogue wave pattern since the system is lin¬ 
ear. However varying the tunneling parameter J changes the 
time at which the peaks appear. If the peaks occur at f = t 
when J - 1, then they would occur at f = t/Jq when J - Jq. 
This scaling can be easily seen from the equation (1). 

Let us now study the underlying mechanism of the discrete 
rogue wave formation in the linear domain. The key idea here 
is the Josephson effect Ell. Our initial wave function (2) is 
uniform but has a position dependent phase. This leads to a 
position dependent current in the lattice. The current in turn 
leads to population imbalance between adjacent lattice sites. 
If the position dependent phase is chosen appropriately, these 
currents add up altogether to form a high density peak at some 
points (the focal points). The Josephson current can easily be 
understood for the exactly solvable two-site system. We follow 
Feynman’s simple two-state phenomenological approach 
to describe Josephson effect. Assume that coupling between 
the the two lattice sites equals J and 'Ll and T '2 represent wave 
functions on the first and second sites, respectively. We fur¬ 
ther assume that the densities at each sides are initially equal. 
Therefore the initial wave functions are given by 'Ll = ^/Noe'^' 
and 'T 2 = where 61 and 62 are the phases on the two 

sides and Nq is the density of particles. This phase difference 
between the lattice sites induces a current, known as Josephson 
current, which is proportional to Jno sin (02 - Oi) ED- Our sys¬ 
tem is indeed a many lattice site generalization of this two site 
system. The Josephson current is position dependent since the 
relative phase of two sites is position dependent. The induced 
Josephson current is symmetrical with respect to the center and 
maximum at the edges of the lattice when C 2 + 0. So, the pop¬ 
ulation imbalance occurs first at the edges and particles start to 
flow towards the center. Under some appropriate conditions, the 
currents at each sites don’t cancel each other and particles are 
focused to some points, where peaks occur. As a result, we say 
that modulational instability is not responsible for rogue waves 
in our system. Instead an analogue of a Josephson-like effect 
can be used to explain the formation of linear discrete rogue 
waves. 

Let us now discuss robustness of rogue waves in the linear do¬ 
main. Suppose the same parameters already used in the Fig-La. 
We consider robustness against noise on the initial density. We 
expect that rogue waves are robust against noise. This is be¬ 
cause Josephson current between adjacent sites is not destroyed 
in the presence of such noise. The Josephson current exists 
even when the system is not uniformly distributed at f = 0. In 
our numerical computation, we add a large amount of noise: 

= A(1 -H 0.3W„)e'®", where VT„ are random numbers with 
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Figure 1: The absolute of the wave function for N = 200, g = 0, A = 0.4 
ni = N/2 and «2 = 0. The system has a single peak when ci = ;r and C 2 = 0.005 
(a) W-shaped structure when ci = 0 and C 2 = 0.03 (b) and zigzag-shaped 
structure when ci = 0 and C 2 = 0.06 ( c). 

zero-mean distribution in the interval [-1,1]. The Fig-2.b dis¬ 
plays time evolution of this initial wave packet. We see some 
fluctuations on the background and a single peak as we ex¬ 
pect. Let us now study robustness against diagonal disorder 
El. We expect a similar picture since random on-site potential 
doesn’t destroy Josephson current, either. As can be seen from 
the Fig.2.a, the pattern formation remains almost the same even 
when the diagonal disorder is large. 

2.7. The effect of Noninear Interaction 

Let us now consider g + Q. Positive and negative values of 
g correspond to an repulsion and attraction. It is well known 
that rogue waves can be formed on a finite background if the 
system is modulationally unstable. In this subsection, our aim 
is not to study rogue wave formation due to modulational insta¬ 
bility. Instead, we analyze the effect of nonlinear interaction on 
the process already present in the linear system. 



Figure 2: The same parameters as in the Fig.l.a are used. The peak is robust 
against even high amount of noise as large as 5 percentage diagonal disorder 
(a) and 30 percentage noise (b). 


During the formation of the peak, the nonlinear interaction 
plays an important role. The nonlinear interaction leads to 
either self-focusing or self-defocusing. In the limiting case 
C 2 = 0, both self-focusing and self-defocusing occur depend¬ 
ing on the values of ci at fixed nonlinear interaction. This was 
experimentally achieved in Il29l . Analogously, we also expect 
self-focusing and self-defocusing even when C 2 0. Let us 
now introduce attractive nonlinear interaction to the system de¬ 
scribed in the Fig-La. If C 2 = 0, we would see self-defocusing 
since ci - n Il29l . We also numerical see self-defocusing when 
C 2 = 0.005. The peak density is lowered by the attractive inter¬ 
action and the width consequently increases. If attractive inter¬ 
action is strong enough, then the peak is completely destroyed 
as can be seen in the Fig.3.a. In other words, rogue wave is not 
formed for strong attractive nonlinear interaction. We would 
like to note that rogue waves given in the Fig.Lb and Fig.l.c. 
are destroyed at much stronger attractve nonlinear interaction. 
The physics of rogue wave is drastically changed if repulsive 
nonlinear interaction is introduced. If g is big enough, then the 
peak doesn’t vanish because of the self-trapping mechanism as 
can be seen in the Fig.3.b. Chaotic fluctuations occur on the 
background while the rogue wave survives with a slight change 
on its amplitude in time. The discrete rogue wave that appears 
from nowhere and survives for a long time was studied in 1^ . 
Note that self-trapping depends on c\ and C 2 . For example, we 
don’t numerically see long-living discrete rogue waves in the 
presence of repulsive nonlinear interaction if the same parame¬ 
ters in the Fig.3.b and Fig.3.c are used. 


3 





Figure 3: The same parameters as in the Fig. La but g = -1 (a) and g = 1 (b) 
are used. The single peak is destroyed when the attractive nonlinear interaction 
is strong enough and long-living discrete rogue wave is formed when repulsive 
interaction is strong enough. 


3. Linear Rogue Waves 



Figure 4: Linear continuous rogue wave for a system described by the 
Schrodinger-like equation. The parameters are given hy n = 0, co = u= I, 
and Ai = A 2 = 1. 

amplitude of the rogue wave can be enhanced by varying A 2 at 
fixed Ai. Note that the amplitude of the wave function becomes 
maximum when L(t) is minimum. This wave can be considered 
as a linear rogue wave since it’s amplitude can be much larger 
than the background amplitude. Furthermore it appears from 
nowhere and vanishes without a trace since = +oo)| — > |Ai|. 
To this end, we would like to write down the most general rogue 
wave solution for the linear Schrodinger-like equation. Since 
the system is linear, a combination of i/'m for different values of 
m is still a rogue wave solution. T* = ^ where are 

m 

arbitrary constants. This is the most general form of the rogue 
wave solution for the Schrodinger-like equation. 


So far, we have analyzed rogue wave formation in a linear 
discrete lattice. We have discussed that Josephson effect plays 
a vital role on the formation of discrete rogue waves. Now, a 
question arises. Is there a linear rogue wave in a continuous 
system described by Schrodinger-like equation? In this section, 
we will answer this question. Intuitively, we expect such waves 
since a continuous system can be thought of an infinitesimally 
spaced discrete system with infinitely many sites. Let us con¬ 
sider the Schrodinger-like equation 


.5iA _ 
dt dx^ 


(3) 


The system admits an exact analytical rogue wave solution 


if/m = exp {-iFn - 

VL 


wx"., a/wx 


2L 


)(4) 


where the quantum number m - 0,1,2,... is a positive inte¬ 
ger, AI and A 2 are constant amplitudes, k is the wave vector, 

E - k , F„ — (n + 1/2)—^ and are the Hermite poly- 
_2 L^ _ 

nomials, L{t) = + w^/4 - 2ut + fi, u and w are arbitrary 

constants and dot denotes time derivation. This solution can 
be verified just by direct substitution of this solution into the 
Schrodinger-like equation. 

Let us analyze this exact solution. In the Fig.4, we plot the so¬ 
lution when m — 0. Strikingly, we see a doubly localized peak 
with a high amplitude. We would like to emphasize that the 


4. Conclusion 

In this paper, we have predicted the existence of both contin¬ 
uous and discrete linear rogue waves. The underlying physics 
for the formation of linear rogue waves can not be explained 
by modulational instability. We have discussed that Joseph¬ 
son effect can be used to explain linear discrete rogue waves. 
We have presented the exact linear rogue wave solution for the 
Schrodinger-like equation. It is tempting to find an exact rogue 
wave solution (not triggered by modulational instability) for the 
nonlinear Schrodinger equation. Finally, we say that an experi¬ 
mental realization of our predictions can be achieved with cur¬ 
rent technology. 
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